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Abstract 

Context. Theoretical and numerical studies of star formation have shown that a magnetic field can greatly influence both disk forma- 
tion and its fragmentation, with even relatively low magnetic field strengths being able to prevent these processes. However, very few 
studies have investigated the combined effects of a magnetic field and turbulence. 

Aims. We study the collapse of turbulent, magnetized prestellar cores, focusing on the efi"ects of magnetic diffusion, and misalignment 
between rotation axis and magnetic field, on the formation of disks, fragmentation and the generation of outflows. 
Methods. We perfom three-dimensional, adaptive mesh, numerical simulations of magnetically super-critical collapsing dense cores 
of 5 Mq using the magneto-hydrodynamic code Ramses. A turbulent velocity field is imposed as initial conditions, characterised by a 
Kolmogorov power spectrum. Different levels of turbulence (a laminar case as well as subsonic and supersonic cases) and magnetiza- 
tion (from weak to strong magnetization) are investigated, as well as three realisations for the turbulent velocity field. 
Results. The turbulent velocity field imposed as initial conditions contains a non-zero angular momentum, which is responsible for 
a misalignment of the rotation axis with respect to the initial magnetic field, and an effective turbulent diffusivity in the vicinity of 
the core. Both effects are responsible for a significant decrease of the magnetic braking, and facilitate the formation of early massive 
disks. These disks can fragment even with /j ~ 5 at late times, in contrast to simulations of 1 Mq cores, where fragmentation is 
prevented for such values of fi. Slow asymmetric outflows are always launched, and they carry a mass which is comparable to that of 
the adiabatic first core. 

Conclusions. Because of turbulence-induced misalignment and magnetic diffusivity, massive disk formation is possible; nevertheless 
their mass and size are much reduced compared to those of disks formed in unmagnetized collapsing cores. We find that for jj. > 5 
fragmentation can occur 

Key words. Magnetohydrodynamics (MHD) - Stars: low mass, formation 



1. Introduction 

Disks are a key element in star formation: protostars proba- 
bly grow by accreting material from protostellar accretion disks 
(iLarson 20031), which at later stages are the natural progeni- 
tors of planets (Lissauer 1993). The formation of disks is a 
challenging theoretical question, and the conditions for their 
presence remain unclear. Magnetic fields certainly play an im- 
portant role, and can efficiently transport angular momentum 
away from the central regions of the collapsing core, prevent- 
ing the early formation of massive disks. This occurs even 
for relatively low rnagnetic field intensities (iMellon & Lill2008t 



iPrice & Batel 120071; iHennebelle & Fromand |2008|) . at a level 



largely in accordance with observational measurement of core 
magnetizations(Crutcher 1999; Falgarone et al. 2008). However, 
this so-called magnetic braking "catastrophe" may in fact be 
compatible with observations. Unlike Class I and more evolved 
sources, the presence of disks around Class sources is not well 
established, and recent studies show no clear evidence of disks 
or fragmentation in Class protostars (Maurv et al. 2010). 

Several studies have investigated non-ideal magnetohydro- 
dynamics (MHD) effects on the magnetic braking and on disk 
formation. While Mellon & Li (2009); Li et al. (2011) found 
that ambipolar diffusion does not allow the building of massive 
disks, the effects Ohmic dissipation remain less clear, with only 



small disks able to form(~ 10 AU iDapp & Bas u"2010^. Other 
means of redistribu t ing m agnetic flux were also investigated. 
iKrasnopolskv et al] (1201 2|) found that flux expulsion from the 
inner core regions by both ambipolar and Ohmic diffusion, can 
promote the development of an interchange which is very effec- 
tive in advecting outward the magnetic field, nevertheless they 
did not observe disks. 

The effects of turbulence were investigated by 

ISantos-Lima et al ] (120 12h who found that an effective tur- 
bulent diffusivity (of the same order of mag nitude as the 
enhanced resistivity of I Krasnopolskv et al.ll201Cil) is sufficient to 
decreases the magnetic braking efficiency, and g enerate disks. 
Simila rly, in the context of massive star formation. ISeifried et al.l 
(1201 2h argued the lack of coherent rotational velocity field in 
their turbulent simulations prevents the build up of a toroidal 
magnetic field component, and reduces the magnetic braking. 
However, they also noted that turbulence induces a strong mis- 
alignment between the rotation axis of the disk and the magnetic 
field. This misalignment was shown in fact to be critical in the 
trans port of ang ular momentum by Hennebelle & Ciardi (2009h 
and Joosetal. (2012). In ideal MHD, Henneb elle & CiarS 
(2009) and Joos et al. (2012) investigated the influence of the 
orientation of the rotation axis with respect to the magnetic field 
on disk formation, showing both analytically and numerically 



1 



Marc Joos et al.: Influence of turbulence on low-mass star formation 



that the tilt of the rotation axis can reduce significantly the 
magnetic braking, allowing the formation of massive disks. 

Understanding the formation of multiple systems is another 
crucial issue related to the angular momentum distribution. A 
large fraction of stars are actually observed to be in binary or 
higher multiple systems (Duquennov & Mavor 1991). The ques- 
tion of the multiple system formation was addressed by several 
authors using numerical simulations. While the magnetic field is 
a key ingredient of star formation theory, simulations investigat- 
ing multiple star formation were until recently mainly hydrody- 
namical. Fragmentation of collapsing cores is thought to be the 
most plausible process to form m ultiple systems. Pioneer work 
from lBoss & Bodenheimerl (11979 ) already suggested that bina- 
ries could form from fragmentation of the collapsing cloud. The 
robustness of these results, as well as the initial parameters of 
the collapse (as the initial shap e of the cloud) , were thoroughly 
investigated (see for example iBonnell et al.l Il99ll: iBate et all 
11995). Some authors took into account the effect of turbulence 
(Goodwin et al. 2004a b), showing that binary and multiple sys- 
tem can form even with low levels of turbulence, and without a 
large initial amount of rotational e nergy. Fragmentation of the 
protostellar disk was also studied (IStamatellos et al.ll2007h . il- 
lustrating that companion stars could form in an extended mas- 
sive disk around the protostar. The eff ects o f radiative transfer, 
without OOffner et all 12009) or with (lOffneret all 120 10) turbu- 
lence were also addressed. The main conclusion being that while 
radiative transfer inhibits fragmentation in the disk, turbulence 
tends to favour fragmentation of the cloud. 

Few studies have investigated the influence of a mag- 
netic field on the fragmentation process, and they reveal in 
particular that it can very efficiently prevent fragmentation 
(Hennebelle & Teyssier 2008; Machidaetal. 2005). The mag- 
netic pressure has a stabilizing influence, particularly because 
of the growth of the toroidal component of the magnetic field 
which is responsible for disk stabilization. Radiative processes 
do n ot enhance fr agmentation in the presence of a magnetic 
field (| Comme rcon et al. 2010), on the contrary, the more realistic 
heating of the disk inhibits the fragmentation even further. Fewer 
studies investigated the role of turbulen ce in magnetized low- 
mass cores, as for example in the work of iMatsumoto & Hanawal 
( 1201 1 ), whose simulations, although not really focused on frag- 
mentation itself, showed no signs of it. 

The goal of this paper is to study in detail disk formation and 
fragmentation in the presence of both a magnetic field and turbu- 
lence. The effects of turbulence highlighted in the present work 
are twofold. First, turbulence generates an effect ive magnetic 
diffusivity by promoting magnetic reconnection (IWeissI [l966t 
[Santos-Lima et al. 2012), which helps to expel the magnetic flux 
from the central regions of the core. Second, it is responsible for 
a "local" misalignment between the rotation axis and the mag- 
netic field lines, which results in a signifi cant decrease of the 
magnetic braking efficiency, as discussed in I Joos et"aL l (l20Tl . 

The plan of the paper is the following. In section 2, we 
present our numerical setup and initial conditions. The effects 
of turbulence are presented in section 3: the diffusion of the 
magnetic field, the impact on the orientation of the angular mo- 
mentum with respect to the magnetic field, and its implications 
on the magnetic braking efficiency. The consequences on star 
formation are presented in section 4: the effects on disk forma- 
tion, fragmentation, and outflows. In section 5, we discuss the 
influence of misalignment on a laminar case and present results 
for different realisations of the turbulent velocity field. Section 6 
concludes the paper. 



2. Numerical setup and initial conditions 

2.1. Numerical setup 

We perf orm 3D numerical simulations wit h the AMR code 
Ramses (l Tevssieiil2002l: iFromang et al.ll200 6'). Ramses can treat 
ideal MHD problems with self-gravity and cooling. The mag- 
netic field evolves using the constrained transport method, pre- 
serving the nullity of the divergence of the magnetic field. The 
high resolution needed to investigate the problem is provided 
by the AMR scheme. Our simulations are performed using the 
HLLD solver (Miyoshi & Kusano 200l). 

The calculations start with 128^ grid cells. As the collapse 
proceeds, new cells are introduced to ensure the resolution of the 
Jeans length with at least ten cells. Altogether we typically use 
10 AMR levels during the calculation, providing a finest spatial 
resolution of ~ 0.4 AU. 

2.2. Initial conditions 

We consider simple initial conditions consisting of a spherical 
cloud with a mass of 5 Mq. Compared to the one solar mass 
clouds usually studied in the context of low-mass star for mation, 
more massive clouds are expected to be more turbulent ( LarsonI 
198jJ. The density profile of the initial prestellar core is given 
by 



with po - 1.2 X 10" g cm" the central density, kq - 2.3 x 
10"^ pc the radius of the central region of the prestellar core, 
where the density is approximately uniform, and = 0.14 pc 
the boundary of the prestellar core; these are in accordance with 
observations (Andre et al. 2000; Belloche et al. 2002). The ratio 
of the thermal over gravitational energy is about 0.25 whereas 
the ratio of the rotational over gravitational energy /? is about 
0.03. The magnetic field intensity is quantified through the mag- 
netization parameter fi, defined as the mass-to-flux over the crit- 
ical mass-to-flux ratio (its critical value being ^ l/(27rG'^^)). 
Various magnetizations are studied, corresponding to a magne- 
tization parameter fi - 2, 5 (magnetized super-critical cloud, in 
agreement with observations, as pointed out in the introduction) 
and 17 (very super-critical cloud). The magnetic field is initially 
aligned with the rotation axis. We emphasize that the /i = 5 
is an intermediate case: in the /i = 2 case, the magnetic field 
dominates the dynamics whereas the fj. - 17 case is a quasi- 
hydrodynamical case. Therefore in the rest of the paper particu- 
lar attention will be payed to the yu = 5 case. 

A random phase velocity field with a Kolmogorov power 
spectrum is imposed as initial condition. The turbulent velocity 
field is a mixture of compressible and solenoidal modes gener- 
ated, nearly in equipar t ition, using a fractional Brownian mo- 
tion, following iLevrie 3 (12004 . The largest wavelength excited 
is comparable to the size of the central region of the prestel- 
lar core, rs- We explore mostly one realisation of the velocity 
field, but we run two other realisations for the fi = 5 simulations 
in order to assess the robustness of our conclusions. Prestellar 
cores are thought to lay in a transonic turbulent environment 
(IWard-Thompson et al.l 120071) . Therefore, various initial turbu- 
lent energy are studied: eturb = 0, 0.1, 0.2 and 0.5, where etuA 
is the ratio of the turbulent over gravitational energy /Sturb/^^grav 
These correspond to a laminar case, a subsonic case, and two 
supersonic cases respectively. Specifically, the values eturb =0.1 
corresponds to a Mach number Ai - 0.85, eturb = 0.2 to Al = 1 .2, 
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Table 1. List of performed simulations, with /i the magnetization 
parameter, etuA - Emih I £grav the ratio of turbulent over gravita- 
tional energies, and the realisations performed. 
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Figure 1. Specific angular momentum (upper panel) and initial 
angle between the rotation axis and the .T-axis at f = yr. The 
shaded area corresponds to the central region of the prestellar 
core of radius tq. 



and fturb - 0.5 to = 1.9. Table [T]lists all the simulations pa- 
rameters. 

In addition to the angular momentum imposed initially on 
the cloud, the turbulent velocity field can generate a supplemen- 
tary non-zero angular momentum. Figure [T(a)]shows that in the 
inner parts of the prestellar core, the specific angular momen- 
tum is larger when turbulence is present. In a sphere of radius 
ro, the angular momentum is larger by a factor ~1.3 to ~2, and 
the larger the turbulent energy, the larger the angular momen- 
tum. However outside the central region of the prestellar core, 
the angular momentum tends to converge to a unique value, in- 
dependently of the initial turbulent energy. This also implies the 
presence of an initial tilt of the rotation axis with respect to the 
jc-axis, which is given by the initial direction of the magnetic 
field, as shown by Fig. |l(b)| With turbulence, this angle is ini- 
tially larger than 30°. However, due to the turbulent motions, 
this angle evolves rapidly after the beginning of the simulation, 
as discussed in the next section. 

To avoid the formation of a singularity and to mimic the fact 
that at high density the gas becomes opaque i.e. nearly adiabatic, 
we use a barotropic equation of state 



1+1^ 

Pad 



2/3 



(2) 



where p^d is the critical density over which the gas becomes adi- 
abatic; we take pad = lO"'-' g cm"-'. When p > pad the adiabatic 
index y is therefore equal to 5/3, which coiTesponds to an adia- 
batic mono-atomic gas. At lower density, when the gas is isother- 
mal, Pjp is constant, with c^ o ~ 0.2 km s"'. The corresponding 



free-fall time is ffj ~ 60 kyr (for an initial central density of about 
1.2 X 10-^** g cm-3). 

Note that for a 5 Mq prestellar core, it is likely that 
the barotropic equation of state is not accurate enough, 
and th at radiative transfer calculat io ns should be em - 
ployed ("Commerc on et al.ll2010l l2011allbt iTomida et al.ll2010h . 
However, such calculations are computational intensive and do 
not allow an easy exploration of the space parameters. 



3. Effects of turbulence 

In this section, we will discuss the main two main effects of tur- 
bulence on the collapse of a prestellar core: namely that the mag- 
netic field can be efficiently diffused out of the inner regions of 
the collapsing core owing to an effective turbulent diffusivity, 
and that the angular momentum induced by the initial turbulent 
velocity field is responsible for a "local" misalignment between 
the rotation axis and the magnetic field. This two effects have 
important consequence on the magnetic braking, and therefore 
on the angular momentum distribution. 



3. 1 . Magnetic diffusion 

Let us begin by investigating the diffusion of the magnetic field 
in the regions around the first Larson's core (the adiabatic core, 
which will be refeiTed to "first core" in the rest of the paper). 

To estimate the diffusion of the magnetic field in our simu- 
lation, we first estimate the magnetization around the first cores. 
All physical quantities are computed for a sphere of radius R; 
to express the magnetic flux, given by jtR^{B), we consider the 
mean value of the magnetic field on the sphere (B) (which is vol- 
ume averaged). The mass of the sphere is given by M(R) and the 
critical mass-to-flux ratio is defined as (M/<I>)cr = l/(2ny/G). 
The time evolution of the magnetization parameter i-i{t)\R of a 
sphere of a given radius R can therefore be expressed as 



1 



M(R) 



(3) 
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Figure!. Evolution of the magnetization parameter within a radius R - 100, 500 and 1000 AU around the first core (left, central 
and right panels respectively) for an initial yu of 17, 5 and 2 (upper, middle and lower panels) and for etuA = 0, 0.1, 0.2, and 0.5. 
f = yr corresponds to the formation of the first core. 



We compute jj(t)\R for R = 100, 500 and 1000 AU. 

Figure |2] shows the evolution /v(f)l« for initial values of fi - 
17, 5 and 2. Time f = yr corresponds to th e formation of th e 
first core. We note that this is different from I Joos et al.l (l2012h : 
here the choice is motivated by the fact that turbulence modifies 
the gas dynamics and delays the collapse. The actual delay de- 
pending on the realisation of the turbulent velocity field and on 
the turbulent energy. 

For all magnetizations, the magnetic flux is initially well- 
conserved at large scale: within a sphere of 500 or 1000 AU, 
the estimated magnetization parameter remains approximately 
constant. However, at later times, the magnetic flux starts to de- 
crease, for example for R = 500 AU the magnetization parameter 
increases by a factor ~2. 

In the laminar case (eturb = 0), the magnetic flux is first 
well-conserved at smaller radii (within a sphere of radius R - 
100 AU). However at times >5000 yr after the formation of the 
first core, it starts to be diffused up to ~2 to 10 times its initial 
value, this is likely a consequence of the interchange instability, 
as discussed by Kr asnopolsky et al. (2012). 

We find that with an initial turbulent velocity field (eturb - 
0.1, 0.2, and 0.5) the magnetic flux is diffused out more effi- 
ciently in the vicinity of the first core (R < 100 AU). The larger 
the initial turbulent energy, the larger is/i. For example, fi is typ- 
ically ~3 to 10 times larger than its initial value in the subsonic 
case (eturb = 0.1), whereas it is ~10 times larger than its initial 
value for et^-b - 0.2, and ~10 to 30 times larger for eturb = 0.5. 

Following the method of Fromang & Papaloizou (2006), 
we have also computed an approximate estimate of the tur- 
bulent resistivity. We obtain a magnetic diffusivity of about 



2x 10^1 cm^ s-^ for eturb = 0.2, and 4x 10^' cm^ s"' for eturb= 0.5. 
Therefore, the larger the turbulent energy, the larger the diffusiv- 
ity coefficient, which is in good agreement with the measured 
diffusion of the magnetic field. There is roughly a factor 2 be- 
tween the diffusivity coefficient for the eturb - 0.2 and 0.5 cases, 
and consequently the magnetization parameter in the vicinity of 
the first core was increased by a factor two between those two 
turbulent cases. 

We note that even a low level of turbulence is sufficient to ef- 
ficiently diffuse the magnetic flux from the first core. The signifi- 
cant diffusion of magnetic field even at larger radii (R = 500 AU) 
is crucial for the transport of angular momentum, as strong mag- 
netic braking occurs in those regions, and we can expect it to 
influence the formation of protostellar disks. 

3.2. Orientation 

As we discussed in lJoos et al ] (120 12h . the orientation of the ro- 
tation axis with respect to the direction of the magnetic field has 
important consequences on the magnetic braking and the trans- 
port of angular momentum. One of the main conclusions was 
that even with a relatively strong magnetic field (fi > 3), a rela- 
tively small misalignment (> 20°), is sufficient to enable the for- 
mation of massive disks. Although those simulations were per- 
formed without turbulence, we suggested that turbulent motions 
in the ISM cou ld be at th e origin of this misalignment (see the 
section 3. 1 .5 of lJoos et al.l (l2012i) ). In the present case, the turbu- 
lent velocity field imposed as initial conditions induces an initial 
tilt of the angular momentum, with values of the angle laying 
between 20 and 60°. 
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Figure 3. Angle between the angular momentum and the magnetic field for fi - 11 (upper panels), 5 (middle panels) and 2 (lower 
panels), for the four levels of turbulence (eturb = 0, 0.1, 0.2, and 0.5). Their respective orientation is showed in two different regions: 
one defined by 100 < r < 500 AU (left panels) and the other by 500 < r < 1000 AU (right panels). 



Figure [3] shows the time evolution of the angle between the 
angular momentum and the magnetic field for fi = 17, 5, and 
2, and eturb = 0, 0.1, 0.2, and 0.5. The misalignment angle a 
is calculated for two regions: 100 < R < 500 AU, and 500 < 
R < 1000 AU. For smaller radii, the magnetic field is strongly 
twisted around the first core and defining a mean direction of 
the field is unrealistic. While the first region corresponds to the 
disk formation region, we also consider a larger scale since the 
orientation is not as clearly defined in the turbulent cases as in 
the laminar case. 



Figure |3]clearly shows that without turbulence, the magnetic 
field and the angular momentum stay well aligned during the col- 
lapse of the first core and the evolution of the first core. On the 
contrary, with turbulence, the angular momentum and the mag- 
netic field are strongly misaligned. The angle a ranges from 20 
to 80°. In addition, the mean angle changes between the inner 
and the outer regions, indicating that the stracture of the mag- 
netic field is complex. Again, the results for the sub- and super- 
sonic cases are very similar. Even a subsonic turbulent velocity 
field is sufficient to tilt the rotation axis of the first core. 



3.3. Magnetic braking 

To study the transport of angular momentum we estimate the 
magnetic braking in a cylinder defined by a radius R = 300 AU 
and a height H - 300 AU. The vertical and radial components 
of the integrated flux of angular momentum transported by the 
magnetic field are defined as 



FSR) = 
Jo 



B^(r, <f>, +h(r)/2)B,(r, cf>, +h(r)/2) 
An 



and 



FAR) = 

■-In nh(R)/2 



(4) 



Jo J-j 



-KR)I2 



B^(R,cf>,z)B,{R,<l>,z)^^ 

R Rdzacb 

An 



where Bi{R,(p,z) = Bi{r - R,(p,z). We note that F~ is the sum 
of the fluxes through the faces defined by Hjl and —Hjl. In 
the following, we study the specific quantities F, /M and F^/M, 
where M is the mass enclosed in the volume of interest. 

Figure |4] shows the evolution of the vertical (upper panels) 
and radial components (lower panels) of the above-mentioned 
fluxes for yU = 17, 5 and 2, and for eturb = 0, 0. 1, 0.2, and 0.5. 
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Figure 4. Evolution of the vertical and radial components of the flux of specific angular momentum transported by the magnetic 
field in logarithmic scale, for jj - 17, 5 and 2 and eturb = 0, 0. 1, 0.2, and 0.5. f = corresponds to the formation of the first core. 



The magnetic braking is globally less eflicient in the turbu- 
lent cases than in the laminar case. For eturb = 0.1, 0.2 and 0.5, 
both components F. and F,- (upper and lower panels of Fig. |4} 
are weaker by a factor of about 2 compared to the laminar case. 
In the subsonic case (etuit =0.1) the flux is stronger than in the 
supersonic cases (etmb = 0.2 and 0.5), although weaker than in 
flie laminar case. In general, the larger the turbulent energy, the 
weaker the fluxes (both components Fr and Fr). 

We also note that after a transitory phase, these fluxes be- 
come almost constant in time. However, while the generation 
of Br and due to turbulent motions is fast, and F^ and Fr 
grow rapidly, in the laminar case their increase is relatively slow 
(-5000 years). 

To summarize, the fluxes of angular momentum transported 
by the magnetic field are weaker in the turbulent cases because 
(i) of the turbulence-induced magnetic diffusivity which reduces 
the magnetic flux in the central regions, (ii) of the misalignment 
between the rotation axis and magnetic field, which decrease the 
magnetic braking efficiency. 

4. Consequences for star formation 

4.1. Disk formation 
4.1 .1 . Disk definition 

We have seen in the previous section that the presence of a turbu- 
lent velocity field can considerably reduce the magnetic braking 
in the inner regions of the first core. We now study its conse- 
q uences on disk for mation by following the analysis developed 
in lJoos et al. which is first briefly recalled. 

We work in the frame of the disk, which is a cylindrical 
frame whose axis is parallel to the direction of the angular mo- 
mentum of the first core. Since disks are expected to be rea- 
sonably axisymmetric, they are defined on concentric and su- 
perposed rings in which density, v elocity, p ressure and magnetic 
field are averaged. As discussed in lJoos et a l. (2012) a single ro- 
tation criterion is not enough to define disks. We use instead a 
combination of five criteria, which are: 

1. the first criterion ensures that rotational support is signifi- 
cant; > fVr', 

2. the second criterion ensures that the disk is reasonably close 
to hydrostatic equilibrium: > fv{. 



3. to check whether rotational support dominates thermal sup- 
port, we require: py^ > fPth', 

4. a connectivity criterion; 

5. a density criterion is used to avoid large spkal arms: p > 

\(f cm-\ 

For the first three criteria, we take a value f -2. 
4.1 .2. Mass of tine disk 

The mass of the disk as a function of time, for different magne- 
tizations and turbulence levels, is presented in Fig. |5] For each 
case f = corresponds to the time when the first core forms. 
Because of the weak magnetic braking, disks always form for 
fi - n. A result that is not changed significantly by turbulence. 
The observed rapid decrease in the masses of the disks is caused 
by fragmentation, which is discussed in the next section. For 
fi = 5 massive disks can form in the presence of turbulence, their 
masses growing up to 1.0 Mq, which represents more than 20% 
of the prestellar core's mass. This is in contrast with the laminar 
(and aligned) case, where only relatively smaller disks are able 
to form. For the more magnetized case, - 2, massive disks do 
not form at all, with or without turbulence, instead disks with 
masses in the range of about 0.3 to 0.4 Mq can form. We note 
that because of numerical convergence issues, these masses are 
probably overestimated by a factor of at least 2 (see appendixlAli. 

In the laminar case, for the intermediate magnetization (/i = 
5), there are two regimes: before ~5000 yr the mass of the disk is 
negligible, and after ~5000 yr it starts to grow with a slope com- 
parable to the turbulent cases (Fig. |5(b)) . We beheve that it is due 
to a symmetry breaking occur ring at this stage, due to the inter- 
change instability discussed in iKrasnopolskv et alJ (l2012l) . This 
symmetry breaking generates turbulent motions in the vicinity 
of the first core and causes a transport of magnetic flux that can 
be seen in the first panel of Fig. |2(b)| where two regimes can 
also be seen before and after ~3000 yr. This magnetic flux diffu- 
sion induces a decrease of the magnetic braking efficiency; more 
angular momentum is then available to build disk. This effect 
probably depends on the non-ideal MHD processes taken into 
account (and on the numerical methods used when these pro- 
cesses are not explicitly treated), and would require further in- 
vestigations. This is consistent with our previous claim that the 
aligned configuration is peculiar. 
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Figure 5. Mass of the disk as a function of time for ji - 11 (Fij 
^.ti^ 5 (Fig.|5(b)]) and^ = 2 (Fig. [5(3)1). 



Figure 6. Mass of the fragments for fj. - 5. The central, most 
massive first core is denoted by the red diamonds and Hne. 



These resuhs show that due to the reduced magnetic braking, 
in the low and intermediate magnetization cases there is a trend 
to form bigger disks when etuit ^ 0. At the same time, it is still 
clear that for increasing magnetic field strength, and thus mag- 
netic braking, disks tend to be smaller More precisely, for ju -2, 
disks with masses greater than 10% of the mass of the prestellar 
core can hardly form, at the time simulated. Importantly, in the 
intermediate magnetization case fi - 5, the mass of the disks in 
the subsonic case (etuit - 0.1) tends to be only slightly smaller 
than in the supersonic cases. Nevertheless the mass of the disks is 
significantly larger than that of the disks built in a non-turbulent 
environment. Therefore, even a low level of turbulence is sufli- 
cient to form massive disks when /i > 5. 

These results depend also on the realisation of the turbu- 
lent velocity field, as discussed in section l572l and show a large 
variance with the realisations. However, the main trends do not 
change. In most of the cases studied, turbulence favours the for- 
mation of more massive disks. 



4.2. Fragmentation 

Massive disks are prone to fragmenting, and to identify and 
study the fragments we use a friend-of-friend algorithm. In this 
method, a fragment is defined by starting from a pre-defined 
over-density, and then finding all cells above a certain den- 
sity threshold which are spatially connected. Contrary to the 
sink particles technique, it is harder with this method to accu- 
rately follow the fragments, and they can also potentially merge. 
However this technique gives sufficiently accurate information 
on the fragmentation process occurring in our simulations. We 
note that to ensure the validity of those results a treatment of the 
second collapse, as well as radiative transfer calculations, would 

be needed. 

As observ e d pre viously by iHennebelle & Tevssied (l2008l) 
and lJoos et al.l (l2012h for 1 M© prestellar cores, we find that for 
ju = 2 no fragmentation occurs. The disk is too small and remains 
stable. Moreover, the magnetic field, and particularly its toroidal 
component, has a stabilizing effect that efficiently prevents frag- 
mentation. On the other hand, the initial turbulent velocity field 
is too weak to provide density perturbations of sufficient am- 
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Figure 9. Mass of the fragments for jj - 17. 
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plitud e that could trigger fragmentation (iHennebelle & Tevssied 
I2008h . 

For fj. - 5, fragmentation takes place when the disks are 
massive enough. This occurs at late time (approximately 20 
000 yr after the formation of the first core) when the mass of 
the disks reaches about 0.7 Mq. Figure |6] shows, for etmb - 01, 



0.2, and 0.5, the evolution of the mass of the fragments: the red 
diamonds (and the red line) correspond to the first core itself 
whereas the blue dots correspond to the disk fragments. One 
dot corresponds to one fragment at a given time. Fragments can 
merge and fragment again, and up to four fragments of about 
IO^^Mq form. The most massive coiTespond to the central first 
core which forms initially with a mass of about few tenths of a 
solar mass. We note that fragmentation occurs for these 5 Mq 
prestellar cores whereas it was not the case for the 1 Mq prestel- 
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lar cores studied in iJoos et al ] (I2012h . Although the number of 
Jean's masses is initially the same, as the collapse occurs, gas 
becomes denser and temperature becomes higher in the 1 Mq 
prestellar cores than in the 5 Mq prestellar cores, simply because 
of the barotropic equation of state. Therefore, 5 prestellar 
cores have more Jeans' masses than 1 Mq prestellar cores. 

As shown by Fig. Q and [8] the fragmentation of disks also 
occurs for ji - 17. The masses of the fragments for = 17 are 
shown in Fig.|9l Typical masses are a few tenth of a solar mass 
for the bigger ones, to lO^^M© for the smaller ones. In the turbu- 
lent cases (eturb ^ 0), the mass of the fragments can be compara- 
ble to the mass of the central object, whereas without turbulence, 
their mass is <30 % of the mass of the first core. There are up to 
four fragments formed in the laminar and subsonic cases, up to 
eight fragments for etmb = 0.2 and two fragments for etmb = 0.5. 
For eturb - 0.5 or higher, fragmentation is quenched because of 
the increased turbulent support. 




4.3. Outflows 
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Figure 10. Maximum of the projected radial velocity for ju = 5, 
for e,urb = 0.2. 



Outflows and jets are one of the most important tracers 
of star formation, being associated with all phases of proto- 
stellar evolution. They are thought to be magneto-centrifugally 
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Figure 11. Maximum of the projected radial velocity for fi - 2, 
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launched ('Blandford & Pavne 1982; Pudritz & Norman 19831 
iPelletier & Pudritz. .1992; ' .Ferreira, .1997. .Uchida & Shibatal 
1985h . and observations seem to confirm that they play an 



important rol e in th e efficient transport of angular momen- 
tum ( Bacciotti et alj 1200 2). The early formation of outflows 
during the collapse of prestellar dense cores was investi- 
gated in both 2D and 3D MHD simulations ('Mell on & ll 
2008; Hennebelle & Fromang 2008; Ciardi & Hennebelld|2010t 
Baneriee et al. 2006; Machida et al. 2008). The formatio n of out- 
flows in a turbulent medium was discussed by Hennebelle et al.l 
( 1201 Ih . however they modelled much more massive prestellar 
cores (100 Mq) than treated here. 

We define outflows as the gas with \ ■ Cr > 0.2 km s"', 
with e^ the radial unit vector in spherical coordinates (the cen- 
ter of the spherical frame being the centre of the first core it- 
self). As a tracer of the outflow structure. Figures [10] and [iT] 
display the maximum of the projected radial velocity in the 
i/z-plane. Despite the complexity of the magnetic field, outflows 
are launched in all simulations, even in the strongly turbulent 
cases, and reach spatial extensions of a few thousand AU. In a 
turbulent environment these outflows are clearly asymmetric and 
less collimated than in the laminar case. Furthermore the out- 
flows are "better-shaped" in the more magnetized cases, where 
the field is more ordered. The mean outflow velocity (mass- 
averaged) in the turbulent cases is about 0.5 km s while peak 



9 



Marc Joos et al.: Influence of turbulence on low-mass star formation 





Figure 12. Total mass in the outflows for yu = 17, 5, and 2, with- 
out turbulence, with eturb = 0.1, 0.2, and 0.5. 



velocities are ~ 5 km s In general, outflow velocities in the 
turbulent cases are about one order of magnitude lower than in 
the non-turbulent case. 

Figure [T2] shows the mass of the outflows for different mag- 
netization parameters. Mass ejection rates are roughly constant 
~ 10"^ Mq yr"' for n - 2 and 5 with turbulence, and about a 
few 10"^ to few 10""* M© yr"' for ji - 17. The total mass in the 
outflows can be as large as a few xlO^'M©, indicating that a sig- 
nificant fraction of the mass of the first core is actually expelled. 
Indeed, as shown in Fig. [13] the ratio of the mass ejected to that 
accreted onto the first core (defined as the mass with a density 
higher than 10'° cm"^^) is between 0.3 to 0.7. This ratio tends 
to be larger with decreasing turbulent energy and magnetic field 
strength, but remains clearly of the same order. 

The fact that this ratio stays of the same order of magnitude i s 
in good agreement with the work of Matzner & McKe j (l2000h . 
who showed that the efficiency of star formation (e ~ Mt/Mq, 
with M, the mass of the star and Mej the mass in the outflows) 
is independent of the mass of the star. This has important con- 
sequence on the relation between the core mass function (CMF) 
and the initial mass function (IMF). If the star formation effi- 



Figure 13. Ratio of the total mass in the outflows and the mass 
of the first core for /i = 17, 5, and 2, without turbulence and with 
etu^b =0.1,0.2, and 0.5. 



ciency does not depend on the mass of the core, the CMF will 
naturally evolve into the IMF with a factor ~ e. This is therefore 
a plausible explanation of the evolution of the CMF into the IMF. 
We note however that we do not find the same numerical ratio 
(between ~1 and 3 instead of ~ 1/3 and 1/2 in Matzner & McKei 
200(1). This is probably because we do not quantify exactly the 
same ratio, as we do not know the final mass of the star. We also 
emphasize that after the formation of the protostar, high velocity 
jets would be triggered, transporting more mass, and therefore 
decreasing the computed ratio. 

5. Discussions 

5.1. Turbulent magnetic diffusivity or misalignment? 

To quantify more precisely the influence of an initial misalign- 
ment angle, a, between the rotation axis and the magnetic field 
on the transport angular momentum, and to disentangle its ef- 
fects from those linked to turbulent diff'usivity, we present here 
a series of simulations without turbulence, but with a 0. In 
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Figure 16. Mass of the disk as a function of time. The cases 
shown are for jj - 5, and include turbulence (eturb = 0.2 and 
0.5), no turbulence and aligned (eturb = and a - 0°), and no 
turbulence but an initial misalignment (a = 57° and a = 80°) 
which are shown with thick lines. 



Figure 15. Time evolution of the vertical and radial components 
of the flux of specific angular momentum transported by the 
magnetic field in logarithmic scale. The cases shown are for 
yu = 5, and include turbulence (eturb = 0.2 and 0.5), no turbu- 
lence and aligned (eturb = and a - 0°), and no turbulence but 
an initial misahgnment (a = 57° and a = 80°). 

particular we take ju = 5, and choose a = 57° and a = 80°, 
which correspond to the typical misalignment angles obtained in 
the simulations with turbulence. We stress however that for tur- 
bulent simulations it is difficult to clearly define an orientation 
angle (cf. 13.21 ). In addition, we remind that the results presented 
earUer for the cases without turbulence, denoted eturb - 0, cor- 
responded to an initially aligned core configuration (i.e. a = 0° 
)■ 

Figure [14] shows the evolution of the magnetization param- 
eter yu(f) given by Eq. [3] within a radius R - 50, 100, 500 and 
1000 AU around the central first core. The magnetization pa- 
rameter increases slightly more rapidly in the misaligned case 
than in the eturb = 0, a = 0° case, but remains significantly lower 
than the turbulent cases. These results show that the diffusion of 
the magnetic field is not greatly affected by misalignment in the 
vicinity of the first core, this is expected since it is the turbulence 
that enhances the magnetic diffusivity. 

Figure [15] shows the evolution of and Fr defined by 
Eq. [4] After a few thousand years these fluxes are lower for 
the misaligned cases than in the the eturb - 0, a - 0° case, 
which indicates that the magnetic braking is less efficient wit h 
an initial misalignment, a result also shown in lJoos et"aL 
Compared to the laminar, aligned case, the magnetic braking 
components also grow more rapidly due to the growth of Br and 
B^. It is also clear that F,- and F, are larger for the cases with 
no-turbulence and an initial misalignment than with turbulence. 
This indicates that while the weaker magnetic braking observed 
in the turbulent cases is partly due to the initial misalignment in- 
duced by the turbulent velocity field, this alone is not enough to 
account for the whole effects of turbulence. 

The effects of misalignment on the mass of the disks are 
shown in Figure [16] In first cores with an initial misalignment 
(and no turbulence), the mass of the disks are very comparable 
to those obtained with turbulence (eturb - 0.2 and 0.5). In par- 
ticular, the a - 57°, eturb = case fits very well the simulation 
with eturb = 0.5, with the disk mass growing to 0.3 Mq about 10 
000 yr after the formation of the first core. The a = 80°, eturb = 
is also in good agreement with the simulation with eturb = 0.2, 



the mass growing to 0.5 Mq about 12 000 yr after the forma- 
tion of the first core. Thus, the tilt of the rotation axis with re- 
spect to the magnetic field plays an important role on the for- 
mation of massive disks in the turbulent cores, although as we 
discussed earlier, turbulence induced magnetic diffusion is also 
playing a role as well. We suggest that a similar situation is also 
likely to occu r in the simulations with turbulence carried out by 
'Seifried et al.l (l2012l) . Although they reported strong misalign- 
ment's (between 60 and 180°, see their Fig. 2, bottom panel), 
the formation of massive disks was attributed to the weak build- 
up of a toroidal magnetic field component. However, we find that 
independently of the magnetization and level of turbulence, the 
observed launching of outflows is a clear signature of the pres- 
ence of a significant toroidal magnetic field component (se e also 
iHennebelle & Fromang|l200l ICiardi & Hennebellell2010l) . and 
as demonstrated here, we attribute massive disk formation to the 
combined eff'ects of turbulent diffusion and misalignment. 

5.2. The effects of different realisations of fhe turbulent 
velocity field 

We investigate here the influence of different realisations of the 
turbulent velocity field on the formation of disks and fragmen- 
tation. These realisations are denoted seedl and seed2, whereas 
the first realisation is denoted seedO. We perform the simulations 
for magnetizations fi - 5 and 2, and for emrb = 0.2 and 0.5. Since 
the turbulent velocity field has a non zero angular momentum, 
these different realisations of the velocity field have different an- 
gular momenta as shown in Fig. [17] and [18] We note in particu- 
lar that for the second realisation of the turbulent velocity field 
(Fig.fTTTi. the angular momentum in the sphere of radius ro (cor- 
responding to the radius of the central region of the prestellar 
core where the density is approximately constant), is the same, 
independently of the turbulent energy. Therefore, the angle be- 
tween the rotation axis and the magnetic field is lower outside the 
prestellar core in this case. For the third reahsation of the turbu- 
lent velocity field, the angular momentum and the initial angle 
are comparable to the ones obtained with the first realisation of 
the turbulent velocity field discussed previously. 

The evolution of the magnetization for the three different 
realisations indicates that the diffusion of the magnetic field 
changes with the turbulent initial velocity field. For example, 
for fi - 5, it is initially more important with the first realisa- 
tion; within 100 AU, the magnetization is already twice its ini- 
tial value when the first core forms. For ji -2, the diffusion also 
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Figure 14. Time evolution of the magnetization parameter within a radius R - 100, 500 and 1000 AU around the central first core 
for an initial fi - 5. The cases shown include turbulence (emrb = 0.2 and 0.5), no turbulence and aUgned (etuib - and a - 0°), and 
no turbulence but an initial misalignment. 
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Figure 17. Specific angular momentum (upper panel) and initial 
angle between the rotation axis and the x-axis at f = yr, for 
the second realisation of the turbulent velocity field. The shaded 
area corresponds to the central region of the prestellar core of 
radius tq. 
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Figure 18. Specific angular momentum (upper panel) and initial 
angle between the rotation axis and the x-axis at f = yr, for the 
third realisation of the turbulent velocity field. The shaded area 
corresponds to the central region of the prestellar core of radius 



varies with the initial velocity field: for eturb = 0.2 it is initially 
comparable, independently of the seed, and it is less efficient in 
removing magnetic flux from the first core regions for the first 
realisation with eturb - 0.5. The angle between the rotation axis 
and the magnetic field also depends on the initial turbulent ve- 
locity field, and for the realisations explored it is a > 20°(which 
is the minimum angle required to form massive disks for ^ - 5, 
iJoos et alj|2012l) . 

The effects of different turbulent realisation on the mass of 
disks is presented in Figure [19] For yU = 5 (Fig. |19(a)| i, mas- 
sive disks always form, but their mass can vary by a factor of 
the order of 2. At f ~ 12 000 yr after the formation of the first 
core, this mass is between 0.3 and 0.6 Mq with etuj-^ — 0.2, and 
between 0.2 and 0.5 Mq with etuA - 0.5. This is partly due to 
the inclination angle, and the turbulent magnetic diffusion, since 
the magnetic field is diffused more efficiently with the first real- 



isation than with the two others. This is also coherent with the 
fact that the initial angular momentum was lower in the second 
realisation of the turbulent velocity field than in the other reali- 
sations. For fi - 2, eturb - 0.5, the trend is different than in the 
other cases: the disk is more massive with the second realisation 
and less massive with the first realisation of the turbulent veloc- 
ity field. We believe the explanation is that the protostar forms 
in the central region of the cloud in the case seedl were it can 
accrete more mass, whereas due to both the strong turbulence 
and magnetic support, it forms in outer regions of the cloud in 
the cases seedO and seed2. 

6. Summary and conclusions 

We have carried out simulations of the collapse of turbulent 
prestellar dense cores of 5 Mq with various magnetization pa- 



12 



Marc Joos et al.: Influence of turbulence on low-mass star formation 



5, Eturb = 0.2 




Eturb =0.5 



(a) 



2, Eturb = 0.2 





Eturb =0.5 




(b) 



Figure 19. Mass of the disk as a function of time for - 5 and 2, for three different realisations, with eturb = 0.2 and 0.5 (respectively 
left and right panel). The three realisations are denoted 'seedO' (solid lines), 'seedl' (dashed lines) and 'seed2' (dotted-dashed lines). 



rameters yu, and presented a detailed analysis of the effects of 
turbulence in the star formation process. In particular, the in- 
fluence of an effective turbulent magnetic diffusivity, and of 
the misalignment between the rotation axis and the magnetic 
field, on magnetic braking. In general, the turbulent velocity 
field imposed as initial conditions has important consequences 
on star formation, in particular regarding the formation of mas- 
sive disks, their fragmentation and the launching of bipolar out- 
flows. We also discussed the influence of the relative orientation 
of the rotation and field, in order to compare its eff'ects to those of 
a turbulent velocity field, and, as a first step of a more complete 
statistical analysis, we performed simulations with two other re- 
alisations of the turbulent velocity field. 
Our main conclusions are the following: 

- Magnetic braking is significantly reduced with turbulence 
for two reasons. First, turbulence efficiently diffuses the 
magnetic field out of the inner regions of the core. The turbu- 
lent magnetic diffusivity is responsible for a decrease of the 
magnetic field strength, and therefore of the magnetic brak- 
ing within the core. Second, the turbulent velocity field con- 
tains a non-zero angular momentum which is responsible for 
a misalignment between the rotation axis and th e magnetic 
field. A s discussed here (and more thoroughly in iJoos et all 
(I2012h l this reduce significantly the magnetic braking. To 
quantify the effects of misalignment we compared the tur- 
bulent results with laminar cases having an initial misalign- 
ment. While this clearly plays an important role in the trans- 
port of angular momentum, the formation of disks and their 
fragmentation, it is not sufficient to account for all the effects 
of an initial turbulent velocity field. 

- Owing to the reduced magnetic braking, massive disks can 
form more easily with turbulence than without. We study 
disk formation using the criteria defined in lJoos et"ar 



and we show that disks as massive as 1 Mq can form after 
25 000 yr with turbulence (in the /i = 5 case), whereas the 
mass of the disk formed in the case without turbulence this 
mass is about 0.5 Mq at a comparable time. For ji - 2, the 
mass of the disks is of the order of 0.4 Mq, which is also 
twice the mass of the disks formed in the laminar case. Due 
to numerical convergence issue, these masses could however 
be overestimated, in particular in the more magnetized case 
(see the discussion in the appendix). As more massive disks 
are able to form, they will be more prone to fragmentation. 
Indeed, the more massive disk fragment into a few clumps, 
with masses growing rapidly from ~ lO"'* to 10"' M© and 
then merge with the central object. We stress however that 
radiative feedback is not treated here and could modify the 
picture. 

- Slow outflows are launched in all simulations. In general tur- 
bulence and the tilt of the rotation axis tend to reduce the 
outflows strength compared to the laminar (aligned) case. In 
turbulent simulations the outflows in the higher magnetiza- 
tion cases (/i = 2) can be stronger than those ejected in the 
lower magnetization cases. In general, the mass of the out- 
flows is a fraction of the mass of the protostellar core, and it 
has a weak dependence on the turbulent energy. 

- A total of three different realisations of the initial turbulent 
velocity field were studied for the cases ^- 5 and jj. -2. The 
results are qualitatively similar: turbulence is responsible for 
an effective magnetic diffusion and misalignment, which re- 
duce the magnetic braking and allow more massive disks to 
form. However, the results are quantitatively different, with 
the various initial turbulent velocity fields giving disks with 
masses within a factor~2 of each other. 

- Regarding the magnetic diffusion and the magnetic braking, 
cores with subsonic turbulence are an intermediate case be- 
tween the laminar and supersonic cases. However, the mass 
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Figure A.l. Mass of the disk as a function of time for ji -2, for Figure A.2. Total mass in the outflows as a function of time for 
ftmb - 0.2 and 0.5 (respectively upper and lower panel), with 10, -2, for etuA - 0-2 and 0.5, with 10, 15, and 20 points resolved 



15, and 20 points resolved per Jeans' length. 



per Jeans' length. 



of the disks built in a subsonic environment are comparable 
to those formed in a supersonic environment. Our analysis is 
therefore robust regarding the turbulent environment of low- 
mass forming stars. 
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Appendix A: Convergence 

We run an additional set of simulations to investigate the con- 
vergence of our simulations. We change the Jeans refinement 
strategy to increase the spatial resolution. A cell was previ- 
ously refined if its size exceeded one-tenth of a Jeans' length 
(Cs(n/Gpy^^). We run simulations with 15 cells per Jeans' length 
(HRl) and another one with 20 cells per Jeans' length (HR2). We 
run this convergence run for ji -2, for etmb = 0.2 and 0.5. 

A.1. Disk 

Figure [ATI shows the mass of the disk for /i = 2, for etmb - 0.2 
and 0.5, with 10, 15, and 20 cells resolved per Jeans' length. It 
shows that convergence is not fully achieved in our simulations. 
Indeed, the disk is less massive by a factor up to ~2 in the higher 
resolution runs. 

A.2. Outflows 

Figure lAT2l shows the mass transported in the outflows for fi -2, 
for eturb = 0.2 and 0.5, with 10, 15, and 20 cells resolved per 
Jeans' length. In the low resolution simulation, the outflows 



transport roughly twice the mass transported in the high resolu- 
tions simulations. 

Even if our conclusions are qualitatively unchanged, these 
convergence runs show that our simulations have not converged 
yet, and that we overestimate the mass in the disk or in the out- 
flows by a factor two for the fJ- - 2 simulations. It is partly due 
to the fact that the disks are small on these simulations; it should 
therefore affect less strongly the simulations with a lower mag- 
netization. 
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